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ABSTRACT 

In this paper, we assemble a catalog of 118 strong gravitational lensing systems from 
SLAGS, BELLS, LSD and SL2S surveys and use them to constrain the cosmic equation 
of state. In particular we consider two cases of dark energy phenomenology: XCDM 
model where dark energy is modeled by a fluid with constant w equation of state param¬ 
eter and in Ghevalier - Polarski - Linder (GPL) parametrization where w is allowed to 
evolve with redshift: w{z) = WQ-\-wijj^. We assume spherically symmetric mass distri¬ 
bution in lensing galaxies, but relax the rigid assumption of SIS model in favor to more 
general power-law index 7 , also allowing it to evolve with redshifts 7 ( 2 ). Our results 
for the XCDM cosmology show the agreement with values (concerning both w and 7 
parameters) obtained by other authors. We go further and constrain the GPL param¬ 
eters jointly with 7 ( 2 ;). The resulting confidence regions for the parameters are much 
better than those obtained with a similar method in the past. They are also showing 
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a trend of being complementary to the supernova la data. Our analysis demonstrates 
that strong gravitational lensing systems can be used to probe cosmological parameters 
like the cosmic equation of state for dark energy. Moreover, they have a potential to 
judge whether the cosmic equation of state evolved with time or not. 

Subject headings: gravitational lensing — dark energy — galaxies: fundamental param¬ 
eters 


1. Introduction 


Signific ance of strong gravitational lensing in cosmology has been recognized quite early. Sem¬ 
inal work of Eefsdaj (1964) demonstrating possibilities of independent determination of the Hubble 
constant from time delays between images gave momentum to the development of strong gravita¬ 
tional lensing theory. However, it has only recently been possible to determine t he Hubble constan t 
from lensing time delays with the precision competitive with other techniques (|Suvu et al.l 12010 1. 
Currently, however, goals and performance of cosmology go far beyond “the quest for two numbers” 
{Hq a nd qn ~ the deceleration parameter). With discovery of the accelerating expansion of the Uni¬ 
verse (|Riess et al.l Il998l : iPerlmutter et al.l Il999l i , understanding this phenomenon has become one 
of the most important issues of modern cosmology. Because there is still no fully convincing hints 
from the side of the theory, we are left with a pragmatic approach to model this phenomenon as 
the so called dark energy — hypothetical homogeneous fluid phenomenologically described by the 
barotropic equation of state p = wp. This approach provides a particularly suitable link between 
theory and observations. Namely, it encompasses the case of cos mological constaiit K {w = — 1) and 
that of the scalar field. If the scalar field settled in an attractor (jRatra fc Peebleslll988l ). w can now 
be constant (with the fundamental demand that w < —1/3 in order to get the acceleration) — this 
is the so-called XCDM cosmology. It contains ACDM as its special case and indeed all constraints 
on w obtained up to now within XCDM cosmology gave w very close to —1. Therefore now there 
is no debate whether the Universe is accelerating, but rather whether w coefficient evolved in time 
in quite recent epochs (say at z = 1 — 6). It is a reasonable question because - as lucidly remarked 
by Linder - if the scalar field stands behind accelerating expansion, it should have evolved since 
the only other case of the scalar field we know (the inflaton) clearly evolv ed because the infla¬ 
tion e n ded. A very c onvenient parametrization of w{z) has been proposed bv IChevalier fc Polarski 
(I200l|) ; iLindeiJ (120031 ): w{z) = + which is essentially the first-order (linear) Taylor expan¬ 

sion in the scale factor - the true gravitational degree of freedom in FRW cosmology. We will denote 
this case as CPL thereafter. The problem is that if we use all conventional cosmological probes like 
SNIa, BAO or CMBR acoustic peaks, we are facing the degeneracy (colinearity) between wq and 
wi parameters (stemming from the fact that w(z) overall should be negative). Standard rulers (like 
BAO and CMBR) and standard candles (SNIa) when used to constrain fim and w parameters acted 
complementarily in the sense that their respective confidence regions in plane were almost 

orthogonal. No such complementarity has been proposed so far for wq, wi parameters, although 
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it has been noticed bv iLindeii (j2004l i that strong lensing systems are promising in this respect. It 


is one of the most important reasons motivating our work. Earlier attempts to use strong lensing 
systems for constraining parameters of the cosmological model were based on two approaches. First 
was the statistical one based on comparison between empirical distrib ution of image sep arations in 


samples 

(Oeuri et al. 

2012) respectively. Second approach made use of 

galaxy clusters in the role of 

lenses (I 

^aczvhski fc Gorski 

1081; 

Sereno 

2002; 

Meneffhetti et al. 

2005; 

Gilmore & Nataravan 

2009; 

Julio et al. 

20ldi where a single lens (cluster) tvpicallv generates ca. 100 images of ca. 30-40 


different sources (distant galaxies). 


W e use another m et hod which c a n be traced back to the papers of e.g. iFutamase &: Yoshida 


(j200l|); iBiesiadal (120061 ): iGrillo et al.l (|2008l ). However, similarly as with the Hubble constant, it 
is only quite recently when reasonable catalogs of strong lenses: containing more than 100 lenses, 
with spectroscopic as well as astrometric data, obtained with well defined selection criteria are 
becoming available. It is also only recently when our knowledge about structure and evolution 
of early type galaxies allows us to undertake the assessment of such important factor as mass 
density profile. Recent works (|Biesiada, Piorkowska, &: Maledl2010l : ICao et al.ll2012h provided hrst 
successful applications of the proposed method in the context of dark energy, although small samples 
have not allowed for stringent constraints. This encourages us to improve and develop it further. 
On a new sample of 118 lenses compiled from SLAGS, BELLS, LSD and SL2S surveys we constrain 
the cosmic equation of state in XCDM cosmology and in GPL parametrization. Moreover we relax 
the rigid assumption of SIS model to general power-law density profile allowing the power-law index 
to evolve with redshifts. 

This paper is organized as follows. In Section [2l we briefly describe the methodology. Then, 
in Section [3] we introduce the strong lensing data used in our analysis. The results are presented 
in Section 0] and concluded in Section [5l 


2. Methodology 


As one of the successful predictions of General Relativity in the past decades, strong gravi¬ 
tational lensing has become a very important astrophysical tool allowing us to use individual lensing _ 

galaxies to measure co smological parameters ( Treu et '^1201)6 : Grillo et ahlboDS : Biesiada. Piorkowska. fc Malec 


2011)1 : iGao et al.l 120121 ). For a specific strong lensing system with the intervening galaxy acting as 


a lens, the multiple image separation of the source depends only on angular diameter distances to 
the lens and to the source, as long as one has a reliable model for the mass distribution within the 
lens. 


Moreover, compared with late-type and unknown-type counterparts, early-type galaxies (ETGs, 
ellipticals) are more likely to serve as intervening lenses for the background sources (quasars or 
galaxies). This is because such galaxies contain most of the cosmic stellar mass of the Universe. 
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This property affects statistics of gravitational lensing phenomenon which leads to a sample domi¬ 
nated by early-type galaxies Isee lKochanek et alJ (j2000|) and references therein). Here, we encounter 
a critical issue that still needs to be addressed in a systematic way: such properties of early-type 
galaxies like their formation and evolution are still not fully understood. Fortunately, a sample 
of well-defined strong lensing systems is now available and can be used to study the mass density 
distribution in early-type galaxies. 

In principle, the lens model often fitted to the observed images is base d on a singular isotherma l 
ellipsoid (SIE) model, in which the projected mass distribution is elliptical (jRatnatunga et al.lll999l j. 
In this paper, we will take a simpler approach and assume spherical symmetry. For a moment we 
will refer to the singular isothermal sphere (SIS) model instead which will then be generalized. 


The main idea of our method is that formula for the Einstein radius in a SIS lens 


6 e = dvr 


a 


SIS 


Dis 


D. 


( 1 ) 


depends on the cosmological model through the ratio of (angular-diameter) distances between lens 
and source and between observer and lens. The angular diameter distance in flat Friedmann- 
Robertson-Walker cosmology reads 


D{z;p) = 


1 


1 Z 


f 


dz' 


h{z';p) 


( 2 ) 


where Hq is the present value of the Hubble function and h{z; p) is a dimensionless expansion 
rate dependent on redshift z and cosmological model parameters are: p = for XCDM 

cosmology or p = for the CPL parametrization of evolving equation of state. More 

specifically, h‘^{z] p) = 0^(1 + z)^ + (1 — f2m)(l + for XCDM cosmology and h‘^{z] p) = 

Hm(l + z)^ + (1 — f2m)(l + 2;)3(i+«’o-i-^i’i) gxp (—^(^) for CPL parametrization. 

Provided one has a reliable knowledge about the lensing system, i.e. the Einstein radius Be 
(from image astrometry) and stellar velocity dispersion asis (form central velocity dispersion cjo 
obtained from spectroscopy), one can use it to test the background cosmology. This method is 
independent of the Hubble constant value (which gets canceled in the distance ratio) and is not 
affected by dust absorption or source evolutionary effects. It depends, however, on the reliability 
of lens modeling (e.g. SIS assumption) and measurements of do- Hopefully, starting with the Lens 
Structure and Dynamics (LSD) survey and the more recent SLAGS survey, spectroscopic data for 
central parts of lens galaxies became available allowing to assess their cent ral velocity dispersions . 
Ther e is a growing evi dence for homologous structure of late type galaxies (IKoopmans et al.ll2006l . 
20091 : iTreu et al.ll2006ll supporting reliability of SIS assumption. 


In our method, cosmological model enters not through a distance measure directly, but rather 
through a distance ratio 

Dis{p) L 


dz 


V^^{zuZs]p) = 


)zi h{z'-,p) 


o.(p) /c 


dz' 


(3) 


0 h(z'\p) 
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and respective observable counterpart reads 


'jyobs _ 


c^Oe 

Airal 


From one side this circumstance is a complication, but it also offers an advantage. Namely, 
as we already mentioned, wq and wi parameters are intrinsically anti correlated, which makes 
it difficult to constrain them. In the T) ratio, however a competition between two angular di- 
ameter distances may lead to positive correlations between wq and wi as discussed recently in 


e^ 


m 


( 120121 1 . 


Biesiada. Piorkowska. & Malec 

(20101: 

Biesiada. Malec & Piorkowska 

(2011 

); 

Cao et al. 


One objection one might rise towards the proposed method is the rigid assumption of the 
SIS model for the lens. Although there are some argume nts that inside Einstein radii total mass 
density follows isothermal profile (|KooDmans et ahl 120061 1. one can expect the deviation from the 
isothermal profile and its evolution with redshift. The last point theoretically stems from the 
structure format ion theory (ellipticals as a result of mergers) and to some extent ha s been supported 


observationally (jR.uff et al. 


2011 


Brownstein et al.ll2012l : ISonnenfeld et al.ll2013bl ). 


Therefore we have to generalize the SIS model to spherically symmetric power-law mass dis¬ 
tribution p ~ r~^. First, we recall that location of observed images, hence the knowledge of Oe 

E^cr where: Re = OeDi is 

D, 


provides us with the mass Mi^^ns inside the Einstein radius: Mi^^ns = '^R^e’^cv where 


the physical Einstein radi us (in [kpc]) in the len s plane and Scr = critical projected 

mass density for lensing ([Schneider et al.lll992l l. Finally we have: 


Miens — 


DgDi 2 


(4) 


4(T A 

If one has spectroscopic data providing the velocity dispersion Uap inside the aperture (more pre¬ 
cisely, luminosity averaged line-of-sight velocity dispersion), then after solving spherical Jeans equa¬ 
tion (assuming that stellar and mass distribution follow the same power-law and anisotropy van¬ 
ishes ) one can assess the dynamical mass inside the aperture projected to lens plane (jKoopmans et al 
20051 ) and scale it to the Einstein radius 

2-7 

fil) 


M, 


TT 


dyn 


^O'apRE 


-^alpDiOE 
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ap 

^ap J 
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where 


/(t) = - 


1 (5-27)(1-7) r(7-I) 


vr 


3-7 


r(7 - 3/2) 


r(7/2 - 1/ 2) 
r(7/2) 


n 2 


( 6 ) 
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By combining Eq. Q and Eq. m, we obtain 


Now, our observable is 








( 8 ) 


and its theoretical counterpart (the distance ratio) {zi , ; p) is given by Eq. ([3|) . 


Eractional uncertainty of P (after calculating all relevant partial derivatives and simplifying 
terms) is 

= ^A{5aapY + {l-^Y{59EY (9) 


5V = 


AP 


P 


Following the SLAGS team we took the fractional uncertainty of the Einstein radius at the level of 
5% i.e. 59e = 0.05 - the same for all lenses. 

Theoretically for a single system one could use Uap, but because we deal with a sample of lenses, 
we shall also transform all velocity dispersions measured within an aperture to velocity dispersions 
within circular a pertur e of radius Reff/2 (half the effective radius) following the prescription of 
Jorgensen et al.l (|l995al lbll: uo = aap{9eff/{29ap))~^'^'^- In the literature it has also been denoted 
as cre 2 . We adopt the convention to denote observed angular Einstein radius or effective radius 
as 9e or 0 e// with the notation Re or Reff reserved for physical values (in [kpc]) of respective 
quantities. In our analysis however, we will consider two cases: with aap and (Tq. While using uo, 
the equations Eq. Eq. ([ 8 ]) and Eq. ([9]) should be modified by replacing Uap with cro- In principle 
(i.e. within the power-law model), the use of aap should work because we rescale to the Einstein 
radius anyway. However, the use of (Tq makes our observable p°^^ more homogeneous for the sample 
of lenses located at different redshifts. One may worry about additional uncertainties introduced 
into the error budget this way since the measurement of 0e// bears its own uncertainty. Arguing 
that fractional uncertainty for the effective radius is at the level of 5% (like for the Einstein radius) 
and bearing in mind that in the Jprgensen formula 0e// is raised to a very small power one can 
estimate that uncertainties of the effective radius contribute less than 1 % to the uncertainty of do. 


Then using routines available within CosmoMC package (jLewis fc Bridlell2nn2l ). we preformed 
Monte Carlo simulations of the posterior likelihood £ ~ exp (—y^/ 2 ) where 


= 


118 

E 

i=l 


l^'Robs 


( 10 ) 


In our fits mass density power-law index 7 was taken as a free parame ter fitted together with 
cos mological parame t ers p . It has been sugg e sted b y iRuff et al.l (j201l|) and further supported 


by iBrownstein et al.l (j 2012 l l: ISonnenfeld et al.l (|2013al L that mass density power-law index 7 of 
massive elliptical galaxies evolves with redshift. On a combined sample of lenses from SLAGS, 
SL2S and LSD they fitted the 7 ( 2 ;;) data with the linear relation and obtained 7 ( 2 :;) = 2 . 12 ]']q'q 4 — 
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0.25]'^qJ 2 ^ + 0.17lQg2('Scatier). Therefore we also performed fits assuming the linear relation: 

li^i) = 7 o + 71 - 2 / treating 70 and 71 as free parameters together with cosmological ones. 


3. Data sets 


In order to implement the methodology described in Section O we have made a comprehensive 
compilation of 118 strong lensing systems from four surveys: SLAGS, BELLS, LSD and SL2S. The 
Sloan Lens ACS Survey (SLAGS) and the BOSS emission-line lens survey (BELLS) are spectro¬ 
scopic lens surveys in which candidates are selected respectively from Sloan Digital Sky Survey 
(SDSS-III) data and Baryon Oscillation S pectroscopic Survey (B OSS). BOSS has been initiated by 
upgrading SDSS-I optical spectrographs (jEisenstein et al.ll201llL The idea was to take the spectra 
of early type galaxies and to look for the presence of emission lines at redshift higher than that of 
the target galaxy. Candidates selected this way were followed-up with HST ACS snapshot imaging 
and after image processing: subtraction of the de Vaucouleurs profile of the target galaxy, those 
displaying multiple images and/or Einstein rings have been classified as confirmed lenses. For 
our purpose SLAG S data comprising 57 strong lenses were taken after iBolton et al.l (12008311 and 


Auger et al. ( 20091 ). And the BELLS data containing 25 lenses were taken from 

(I2OI2I L 


Brownstein et al 


Lenses Structure and Dynamics (LSD) survey was a predecessor of SLAGS in the sense that 
combined image and lens velocity dispersion data were used to constrain the structure of lens¬ 
ing galaxies. However, it was different, because lenses were selected optically (as multiple im¬ 
ages of sources with identified lensing galaxies) and then followed-up spectroscopically. There¬ 
fore, in order to comply with SLAGS and BELLS, we took only 5 most reliab l e lense s from LSD: 
CFR S03.1077: HST 14176-752 26: HST 15433-^5352 af ter iTreu fc KoonmansI (120041 1. Q0047-281 
after Koopmans fc Treu ( 2002 1. and MG2016-I-112 after Treu fc Knnpma,uX ( 2002 ). 


At last the Strong Lensing Legacy Survey (SL2S) is a project dedicated to finding galaxy-scale 
lenses in the Canada France Hawaii Telescope Legacy Survey (CFHTLS). The targets are massive 
red galaxies and an automated Ring Finder software is looking for tangentially elongated blue 


features around (jCavazzi et al. II2014I ). If found they are followed -up with HST arid spect roscopy. 


The data for a total of 31 lenses we include here were taken from ISonnenfeld et al.l (j2013al jbl) 


Our sample is presented in Table 1. It contains all relevant data necessary to perform cos¬ 
mological model fits according to our methodology presented in Section [2l The Einstein radius 
derivation method is more or less consistent across different surveys considered here. After sub¬ 
tracting de Vaucouleurs profile of the deflector, the Einstein radii were measured by fitting model 
mass distributions to generate model lensed images and comparing them to the observed images. 
The mass distribution of the main lenses was modeled as singular isothermal ellipsoids (SIE), while 
in several systems, an external shear component was also added to describe the lensing effect of 
nearby groups or clusters. Even though individual uncertainties of this procedure are different 















































depending on the survey and on whether the image was taken from the Earth or from space (HST) 
there is a consensus that in average the relative uncertainty of the Einstein radius is at the level 
of 5%. Apertures for SLAGS and BELLS were taken as 1.5" and 1" respectively according to the 
source papers. Eor LSD and SL2S lenses we have assessed the aperture size from the sizes x, y 
of the slit reported in source papers. The last column gives the velocity dispersion within a half 
effective radius calculated according to Jprgensen formula. Fig. [T] shows the scatter plot of our 
lensing systems. One can see that inclusion of the SL2S lenses resulted in a fair coverage of lenses 
and sources redshifts. 

Fig.[2]presents histograms of lens and source redshift with the contribution of respective surveys 
over-plotted in the same color convention as in Fig. [TJ Different surveys have the following median 
values of the lens redshifts; SLAGS - zi = 0.215, BELLS - zi = 0.517, LSD - zi = 0.81 and 
SL2S - zi = 0.456. The above values refer to lenses used by us. SL2S survey - an ongoing one - is 
particularly promising for the future since it already reached the maximum lens redshift of zi = 0 . 8 . 


4. Results 


Although the lens redshift zi and source redshift Zg both cover a wide range in our sample, 
distance ratio D is still co nfined to a very compact range of valu es, which leads to poor constraining 
power for Qm parameter (|Biesiada. Piorkowska. &: Maledl2010l L Therefore i n this paper, we fix Qm 
at the best-fit value = 0.315 based on the recent Planck observations (jAde et al.ll2014l L This 
disadvantage of our method, i.e. the necessity of taking a prior for the matter density parameter, 
is to a certain degree alleviated by the benefit of being independent of the Hubble constant which 
cancels in the distance ratio. Gonsequently Hq and its uncertainty do not influence the results. 
Performing fits to different cosmological scenarios on the n = 118 strong lensing sample, we obtain 
the results displayed in Table 2. The marginalized probability distribution of each parameter and 
the marginalized 2D confidence contours are presented in Fig. [3]l6l 


We started our analysis with dark energy phenomenon modeled by barotropic fluid having 
constant equation of state coefficient w and we considered two cases of mass density profiles in lenses: 
a non-evolving power-law density profile and an evolving one (denoted in Table 2 as XCDMl and 
XCDM2, respectively). Power-law exponents 7 , 70 and 71 were treated as free parameters to be 
fitted. 


First, we consider fits on our sample taking the aperture velocity dispersion aap as the lens 
parameter. In the XCDM model where w is the only free cosmological parameter and 7 is the 
mass density power-law index parameter, we obtain w = — l.dblggg and 7 = 2.03 ± 0.06. In the 
second case, when the power-law mass density profile was allowed to evolve: ^{zi) = 70 + 71-^7 the 
best-fit values for the parameters are w = — 1.48lg 94 , 70 = 2.06±0.09, and 71 = —0.09±0.16. One 
can see that w coefficient obtained from the strong lensing sample is consistent with the AGDM 
model {w = —1) at la level. Fits on the 7 parameter also reveal compatibility between our sample 
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of lenses and the previou s smaller combined sample from SLAGS, SL2S and LSD (jRuff et alJl201ll : 


Sonnenfeld et al.M2013bl L In an attempt to constrain cosmology with the CPL parametrization 


describing an evolving cosmic equation of state, we first consider the case that both wq, wi and 
7o, 7i are free parameters (denoted in Table 2 as CPLl). By fitting the CPL model to the full 
n = 118 sample with Uap, we get the marginalized la constraints of the parameters tco = —O.lblJ'gQ, 
wi = —6.95^305 and 70 = 2.08 ± 0.09, 71 = —0.09 ± 0.17. By fixing 7 parameters at our best- 
fit values (denoted in Table 2 as CPL2), the best-fit values for the two cosmic equation of state 
parameters are: wo = —O.lOtj'^g and wi = —6.25'^^''^^. We also show the marginalized la and 2a 
contours of the two par ameter in Fig. [ 6 l It can be seen tha t , comparing to t he previous analysis 
with a smaller sample ([Biesiada, Piorkowska, &: Maled I 2 OIOI : ICao et al.l [2012) , fits for wq and wi 
are significantly improved with a larger strong lensing sample. 


Working on the sample with the aperture corrected velocity dispersion ao, we find that 
the dark energy equation of state parameter (w = —1.15l5'2o ^ ~ very 

well with th e respective valu e derived from Planck observations combined with BAO data (w = 
— I.IS^’JqJq) ( Ade et al.ll2014l ). Constraints on the mass density power-law index parameters (70 = 


2.13^n ?o, 7 i = —0.09 zb 0.17) are even more consistent with the previous analysis ( Ruff et al. 


2011 


Sonnenfeld et al.ll2013bl l especially with uq taken as the lens parameter. In the case of evolving 


cosmic equation of state (CPL parametrization), we obtain wq = — l.OOlJ'gg, wi = — l.Sblgyg and 


7o = 2.14_Q'3g, 71 = —0.10 ± 0.18. By fixing 7 parameters at the best-fit values, the fits for wq 
and wi are: wq = —and wi = —l.Gblggl. One can see fairly good agreement between 
our results obtained by using distance ratio method for strong lensing systems and the concor¬ 
dance ACDM model. In order to compare our fits with the results obtained using sup ernovae la, 
likelihood contours obtained with the latest Union2.1 compilation ((Suzuki et al.ll2ni2l ) consisting 
of 580 SN la data points are also plotted in Fig. 6 . For a fair comparison, the matter density 
parameter fim is also fixed at the Planck best-fit value fim = 0.315 and t he systematic err ors of 
observed distance moduli are also considered in the likelihood calculation (jCao &: Zhull2014l L We 
see that la confidence regions from the two data sets overlap very well with each other. This 
means that the results obtained on the sample of strong lenses are consistent with the SNIa fits. 
This consistency at la level is different from earlier re sults obtained with smaller sample of lenses 
( Biesiada. Piorkowska. fc: Malec 2010 : Cao et al. 2012). 


One can clearly see from Fig [ 6 | that principal axes of confidence regio ns obtained with super- _ 

novae and strong lenses are inclined at high e r angles that in p revious studies ([Biesiada. Piorkowska. fc: Malec 


2010 


; IPiesiada. Malec &: Piorkowska 2011 : Cao et al. 2012). This sustains the hope that careful 


choice of the sample in terms of lens and source redshifts, would eventually realize an ultimate 
dream to ha,v e a complementary prob e breaking the degeneracy in (wo,wi) plane (discussed in 
Linder! (j2004l f: [Piorkowska et al.l (j2013[)). This will be a subject of another study. 


Our method based on distance ratio for strong gravitational lensing systei ns may also contrib ute 
to testing the consistency between luminosity and angular diameter distances ( Cao fc ZhJboila b). 
As it is well known, these two observables are related to each other via Etherington reciprocity 
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relation. Earlier discussions of this is sue can be found in iBassett &: Kuna (|2004l l: lUzan et al 


(200J); iHolanda. Lima Sz Ribeirol ( 2010 1 


In order to study the systematics and scatter in our method we performed the diagnostics 
of residuals. Plots of relative residuals — T>^^)/T>°’^^ as a function of zi, Zg, Re/R eff and 

(To for the case of XCDM cosmology, assuming non-evolving mass density profile in lenses, are 
displayed in Fig. [71 This figure is representative of similar diagnostics for the CPL case and 
evolving mass density profile. One can see that there is no correlation between residuals and the 
redshifts of lenses or sources, as well as with the Einstein radius relative to the effective radius. 
However, there is a noticeable anti-correlation (correlation coefficient ca. -0.6) with the velocity 
dispersion. This trend is especially pronounced for lenses with (Tq < 230 km/s. If one excluded 
small velocity dispersion lenses the result would be comparable t o other diagnostics discussed. 
This effect is probably related to recent finings of IShn et al. (j2014l ) who extended SLAGS survey 
into lower mass region and found that elliptical galaxies with smaller observed velocity dispersions 
are more centrally concentrated. In any case, however there is a considerable scatter left (at the 
level of ±50%). This scatter could be attributable to individual properties of lenses, like their 
environment or deviatio n form the spherical symmetry. Fro m theoretical point of view it has been 
known for a long time ([Falco. Gorenstein fc Shapiro Ill985l ) that there exists a transformation of 
mass distribution in the lens that leaves all observables invariant. This is so called “mass-sheet” 
degeneracy — a topic w hich has recently revived in the context of double Einstein Ring lenses 
([Schneider &: Sluse Il2013l l. However, from the physical point of view one of the main ingredients of 
the “mass-sheet” degeneracy comes from the contamination of secondary lenses (clumps of matter 
along or close to the line of sight). Such factors have not been considered in our methodology and 
their proper inclusion is still challenging. 


5. Conclusions 


In conclusion, our analysis demonstrates that strong gravitationally lensed systems can al¬ 
ready now be used to probe cosmological parameters, espec ially the c osmic equation of s t ate fo r 
dark energy. One may say that the approach initiate d in iBiesiadal (|2006l l: iGrillo et al.l (|2008l L 
Biesiada. Piorkowska. &: Maled ([2010|); ICao et al.l (|20I2l l can be further developed. 


However, there are several sources of systematics we do not consider in this paper and which 
remain to be addressed in the future analysis. Let us start with simplified assumptions underlying 
our method. The first one is related to the interpretation of observed velocity dispersions. In 
this paper, we adopted the spherical Jeans equation to connect observed velocity dispersions to the 


masse s and we did it assuming that anisotropy 0 parameter was zero. As shown in lKoopmans et al 


([20061 1 and in more detailed way in iKoopmans et al.l ([20051 1 (the equations generalizing our Eq. [5l|7| 
to arbitrary /3 can be found there) anisotropy parameter is degenerated with the slope 7 . Therefore 
in our approach power-law index should understood as an effective descriptor capturing both the 
density profile and anisotropy of the velocity dispersions. This effect clearly contributes to the 
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scatter in our results. Another issue is the three-dimensional shape of leasing galaxies, the prolate- 
ness/oblateness of leasi ng galaxies can systematically bias the connection between the mass and 
the velocity dispersion (jChad l2003l b This effect also contributes to the scatter at high redshifts 
and might reveal as a systematic effect at low redshifts. All these effects are hard to be rigorously 
accounte d in context of cosmological studi e s like in this paper. Recent results on this issue can be 
found in iBarnabe. Spiniello fc Koopmans I ( 20141 ). 


The other systematic is the influence of the line of sight (fo reground and background) con- 
tamination. T he problem has been recognized long time ago (jKochanek &: Apostolakisl Il988l : 


Bar-KanalIl996l : iKeeton. Kochanek &: SeliakI 119971 ) with a heuristic suggestion that adding an ex- 
ternal shear to an elliptic lens mode l greatly improves the fits of multiple image configurations. 
Wambsganss. Bode fc Ostrikerl (120051 ) addressed the question of secondary lenses on the line of sight 
and concluded that the role of secondary lenses is a strong function of source redshift and can be 
important in 38% of cases for a source at = 3.5. High redshift sources are advantageous from the 
point of view of dark energy studies and at the same time they are challenging from the point of view 
of lin e of sight contamination. One of the most recent studies (j.laroszvhski fc Kostrzewa-Rutkowska 
2 OI 2 I I trying to quantify the influence of the matter along line of sight on strong lensing used the 
technique of simulations of many multiple image conhgurations using a realistic model of light prop¬ 
agation in an inhomogeneous Universe model (b ased on the Mi ll enium simulation). Further progress 
in this direction has recently been achieved by ICollett et al. (j2014l l in a paper accompanied with 
publicly available code Pangloss. They used a simple halo model prescription for reconstructing 
the mass along a line of sight up to intermediate redshifts and calibrated their procedure with 
ray-tracing through the Millenium Simulation. 

Finally, as proposed bv lRusin, Kochanek fc Keeton ( 2003 1. another approach to constrain the 
power-law index 7 of lensing galaxies is to assume a self-similar mass profile and combine different 
strong lens systems to statistically constrain the mass prohle, which could be combined with the 
method considered in this paper to break the degeneracy between 7 and cosmological parameters. 
Inclusion of this combination into cosmological use of strong lensing systems will be a subject of a 
separate paper. 
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Table 1. Compilation of strong lensing systems 


Name 

Zl 


aap [km/s] 

Oe ["] 

survey 

0ap ["] 

Oeff ["] 

(To [km/s] 

J0151+0049 

0.517 

1.364 

219±39 

0.68 

BELLS 

1 

0.89 

226±40 

J0747+5055 

0.438 

0.898 

328±60 

0.75 

BELLS 

1 

1.24 

334±61 

J0747+4448 

0.437 

0.897 

281±52 

0.61 

BELLS 

1 

2.87 

277±51 

J0801+4727 

0.483 

1.518 

98±24 

0.49 

BELLS 

1 

0.57 

103±25 

J0830+5116 

0.53 

1.332 

268±36 

1.14 

BELLS 

1 

1.1 

274±37 

J0944-0147 

0.539 

1.179 

204±34 

0.72 

BELLS 

1 

1.35 

207±35 

J1159-0007 

0.579 

1.346 

165±41 

0.68 

BELLS 

1 

0.99 

170±42 

J1215+0047 

0.642 

1.297 

262±45 

1.37 

BELLS 

1 

1.42 

266±46 

J1221+3806 

0.535 

1.284 

187±48 

0.7 

BELLS 

1 

0.93 

193±49 

J1234-0241 

0.49 

1.016 

122±31 

0.53 

BELLS 

1 

1.61 

123±31 

J1318-0104 

0.659 

1.396 

177±27 

0.68 

BELLS 

1 

1.06 

182±28 

J1337+3620 

0.564 

1.182 

225±35 

1.39 

BELLS 

1 

1.6 

227±35 

J1349+3612 

0.44 

0.893 

178±18 

0.75 

BELLS 

1 

2.03 

178±18 

J1352+3216 

0.463 

1.034 

161±21 

1.82 

BELLS 

1 

1.35 

164±21 

J1522+2910 

0.555 

1.311 

166±27 

0.74 

BELLS 

1 

1.08 

170±28 

J1541+1812 

0.56 

1.113 

174±24 

0.64 

BELLS 

1 

0.59 

183±25 

J1542+1629 

0.352 

1.023 

210±16 

1.04 

BELLS 

1 

1.45 

213±16 

J1545+2748 

0.522 

1.289 

250±37 

1.21 

BELLS 

1 

2.65 

247±37 

J1601+2138 

0.544 

1.446 

207±36 

0.86 

BELLS 

1 

0.63 

217±38 

J1611+1705 

0.477 

1.211 

109±23 

0.58 

BELLS 

1 

1.33 

111±23 

J1631+1854 

0.408 

1.086 

272±14 

1.63 

BELLS 

1 

2.07 

272±14 

J1637+1439 

0.391 

0.874 

208±30 

0.65 

BELLS 

1 

0.89 

215±31 

J2122+0409 

0.626 

1.452 

324±56 

1.58 

BELLS 

1 

1.76 

326±56 

J2125+0411 

0.363 

0.978 

247±17 

1.2 

BELLS 

1 

1.47 

250±17 

J2303+0037 

0.458 

0.936 

274±31 

1.02 

BELLS 

1 

1.35 

278±31 

J0008-0004 

0.44 

1.192 

193±36 

1.16 

SLAGS 

1.5 

1.71 

197±37 

J0029-0055 

0.227 

0.931 

229±18 

0.96 

SLAGS 

1.5 

2.16 

232±18 

J0037-0942 

0.196 

0.632 

279±10 

1.53 

SLAGS 

1.5 

2.19 

283±10 

J0044+0113 

0.12 

0.196 

266±13 

0.79 

SLAGS 

1.5 

2.61 

267±13 

J0109+1500 

0.294 

0.525 

251±19 

0.69 

SLAGS 

1.5 

1.38 

259±20 

J0157-0056 

0.513 

0.924 

295±47 

0.79 

SLAGS 

1.5 

1.06 

308±49 

J0216-0813 

0.332 

0.524 

333±23 

1.16 

SLAGS 

1.5 

2.67 

335±23 

J0252+0039 

0.28 

0.982 

164±12 

1.04 

SLAGS 

1.5 

1.39 

169±12 
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zL 


Fig. 1.— Scatter plot of our sample of 118 strong lensing systems. One can see a fair coverage of 
redshifts in the combined sample. 
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Histogram of zL 
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Fig. 2.— Histogram of lens (upper panel) and source redshifts (lower panel). Different colors show 
contribution of different surveys superimposed on the overall histogram. The color convention is 
the same as for Fig. [TJ 
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Fig. 3.— Joint fits of mass density slope 7 and w coefficient in the XCDM model. Left panel 
shows the results obtained with velocity dispersion within the aperture aap while on the right panel 
corrected velocity dispersion cto was used. = 0.315 is assumed based on the Planck observations 
(|Ade et al.ll2014l L Marginalized probability density functions for 7 and w are also shown. 



Fig. 4.— Joint fits of w coefficient in the XCDM model and mass density slope parameters ( 70 , 
7 i) in evolving slope scenario 'y(z) = 70 + Jizi. Left and right panels display the results obtained 
by using aap and ao respectively. 



































- 20 - 



Fig. 5.— Joint fits of (wo,wi) evolving cosmic equation of state coefficients in the CPL parametriza- 
tion. Leasing galaxies were assumed to have evolving mass density slope j(z) = 70 + 712 ; and ( 70 , 71 ) 
parameters. Left and right panels display the results obtained by using aap and do respectively. 



Fig. 6 .— Joint fits of (wo, wi) evolving cosmic equation of state coefficients in the CPL parametriza- 
tion. Leasing galaxies were assumed to have evolving mass density slope j(z) = 70 + 71 and ( 70 , 71 ) 
parameters were fixed at our best-fit values in XCDM model. Left and right panels display the 
results obtained by using aap and do respectively. 
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Fig. 7.— Relative residuals in our observable: as a function of lens redshift 

zi, source redshift Zg, Einstein radius relative to the effective radius Re/R eff and the aperture 
corrected velocity dispersion (Tq. 
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Table 1—Continued 


Name 

Zl 


dap [km/s] 

Oe ["] 

survey 

Oa, n 

Oeff n 

do [km/s] 

J0330-0020 

0.351 

1.071 

212±21 

1.1 

SLAGS 

1.5 

1.2 

220±22 

J0405-0455 

0.075 

0.81 

160±8 

0.8 

SLAGS 

1.5 

1.36 

165±8 

J0728+3835 

0.206 

0.688 

214±11 

1.25 

SLAGS 

1.5 

1.78 

219±11 

J0737+3216 

0.322 

0.581 

338±17 

1 

SLAGS 

1.5 

2.82 

339±17 

J0808+4706 

0.219 

1.025 

236±11 

1.23 

SLAGS 

1.5 

2.42 

238±11 

J0822+2652 

0.241 

0.594 

259±15 

1.17 

SLAGS 

1.5 

1.82 

264±15 

J0841+3824 

0.116 

0.657 

225±11 

1.41 

SLAGS 

1.5 

4.21 

222±11 

J0903+4116 

0.43 

1.065 

223±27 

1.29 

SLAGS 

1.5 

1.78 

228±28 

J0912+0029 

0.164 

0.324 

326±12 

1.63 

SLAGS 

1.5 

3.87 

323±12 

J0935-0003 

0.348 

0.467 

396±35 

0.87 

SLAGS 

1.5 

4.24 

391±35 

J0936+0913 

0.19 

0.588 

243±12 

1.09 

SLAGS 

1.5 

2.11 

246±12 

J0946+1006 

0.222 

0.608 

263±21 

1.38 

SLAGS 

1.5 

2.35 

266±21 

J0956+5100 

0.24 

0.47 

334±17 

1.33 

SLAGS 

1.5 

2.19 

338±17 

J0959+0410 

0.126 

0.535 

197±13 

0.99 

SLAGS 

1.5 

1.39 

203±13 

J1016+3859 

0.168 

0.439 

247±13 

1.09 

SLAGS 

1.5 

1.46 

254±13 

J1020+1122 

0.282 

0.553 

282±18 

1.2 

SLAGS 

1.5 

1.59 

289±18 

J1023+4230 

0.191 

0.696 

242±15 

1.41 

SLAGS 

1.5 

1.77 

247±15 

JllOO+5329 

0.317 

0.858 

187±23 

1.52 

SLAGS 

1.5 

2.24 

189±23 

J1106+5228 

0.096 

0.407 

262±13 

1.23 

SLAGS 

1.5 

1.68 

268±13 

J1112+0826 

0.273 

0.63 

320±20 

1.49 

SLAGS 

1.5 

1.5 

329±21 

J1134+6027 

0.153 

0.474 

239±12 

1.1 

SLAGS 

1.5 

2.02 

243±12 

J1142+1001 

0.222 

0.504 

221±22 

0.98 

SLAGS 

1.5 

1.91 

225±22 

J1143-0144 

0.106 

0.402 

269±13 

1.68 

SLAGS 

1.5 

4.8 

264±13 

J1153+4612 

0.18 

0.875 

226±15 

1.05 

SLAGS 

1.5 

1.16 

235±16 

J1204+0358 

0.164 

0.631 

267±17 

1.31 

SLAGS 

1.5 

1.47 

275±17 

J1205+4910 

0.215 

0.481 

281±14 

1.22 

SLAGS 

1.5 

2.59 

283±14 

J1213+6708 

0.123 

0.64 

292±15 

1.42 

SLAGS 

1.5 

3.23 

291±15 

J1218+0830 

0.135 

0.717 

219±11 

1.45 

SLAGS 

1.5 

3.18 

218±11 

J1250+0523 

0.232 

0.795 

252±14 

1.13 

SLAGS 

1.5 

1.81 

257±14 

J1251-0208 

0.224 

0.784 

233±23 

0.84 

SLAGS 

1.5 

2.61 

234±23 

J1330-0148 

0.081 

0.712 

185±9 

0.87 

SLAGS 

1.5 

0.89 

194±9 

J1402+6321 

0.205 

0.481 

267±17 

1.35 

SLAGS 

1.5 

2.7 

268±17 

J1403+0006 

0.189 

0.473 

213±17 

0.83 

SLAGS 

1.5 

1.46 

219±17 
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Table 1—Continued 


Name 

zi 

2^5 

aap [km/s] 

Oe ["] 

survey 

Oap n 

Oeff n 

(To [km/s] 

J1416+5136 

0.299 

0.811 

240±25 

1.37 

SLAGS 

1.5 

1.43 

247±26 

J1430+4105 

0.285 

0.575 

322±32 

1.52 

SLAGS 

1.5 

2.55 

324±32 

J1436-0000 

0.285 

0.805 

224±17 

1.12 

SLAGS 

1.5 

2.24 

227±17 

J1451-0239 

0.125 

0.52 

223±14 

1.04 

SLAGS 

1.5 

2.48 

225±14 

J1525+3327 

0.358 

0.717 

264±26 

1.31 

SLAGS 

1.5 

2.9 

264±26 

J1531-0105 

0.16 

0.744 

279±14 

1.71 

SLAGS 

1.5 

2.5 

281±14 

J1538+5817 

0.143 

0.531 

189±12 

1 

SLAGS 

1.5 

1.58 

194±12 

J1621+3931 

0.245 

0.602 

236±20 

1.29 

SLAGS 

1.5 

2.14 

239±20 

J1627-0053 

0.208 

0.524 

290±14 

1.23 

SLAGS 

1.5 

1.98 

295±14 

J1630+4520 

0.248 

0.793 

276±16 

1.78 

SLAGS 

1.5 

1.96 

281±16 

J1636+4707 

0.228 

0.674 

231±15 

1.09 

SLAGS 

1.5 

1.68 

236±15 

J2238-0754 

0.137 

0.713 

198±11 

1.27 

SLAGS 

1.5 

2.33 

200±11 

J2300+0022 

0.228 

0.464 

279±17 

1.24 

SLAGS 

1.5 

1.83 

285±17 

J2303+1422 

0.155 

0.517 

255±16 

1.62 

SLAGS 

1.5 

3.28 

254±16 

J2321-0939 

0.082 

0.532 

249±8 

1.6 

SLAGS 

1.5 

4.11 

246±8 

J2341+0000 

0.186 

0.807 

207±13 

1.44 

SLAGS 

1.5 

3.15 

207±13 

Q0047-2808 

0.485 

3.595 

229±15 

1.34 

LSD 

1.25 

0.82 

239±16 

CFRS03-1077 

0.938 

2.941 

251±19 

1.24 

LSD 

1.25 

1.6 

256±19 

HST14176 

0.81 

3.399 

224±15 

1.41 

LSD 

1.25 

1.06 

232±16 

HST15433 

0.497 

2.092 

116±10 

0.36 

LSD 

1.25 

0.41 

125±11 

MG2016 

1.004 

3.263 

328±32 

1.56 

LSD 

0.65 

0.31 

347±34 

J0212-0555 

0.75 

2.74 

273±22 

1.27 

SL2S 

0.9 

1.22 

277±22 

J0213-0743 

0.717 

3.48 

293±34 

2.39 

SL2S 

1 

1.97 

293±34 

J0214-0405 

0.609 

1.88 

287±47 

1.41 

SL2S 

1 

1.21 

293±48 

J0217-0513 

0.646 

1.847 

239±27 

1.27 

SL2S 

1.5 

0.73 

253±29 

J0219-0829 

0.389 

2.15 

289±23 

1.3 

SL2S 

1 

0.95 

298±24 

J0223-0534 

0.499 

1.44 

288±28 

1.22 

SL2S 

1 

1.31 

293±28 

J0225-0454 

0.238 

1.199 

234±21 

1.76 

SL2S 

1 

2.12 

233±21 

J0226-0420 

0.494 

1.232 

263±24 

1.19 

SL2S 

1 

0.84 

272±25 

J0232-0408 

0.352 

2.34 

281±26 

1.04 

SL2S 

1 

1.14 

287±27 

J0848-0351 

0.682 

1.55 

197±21 

0.85 

SL2S 

0.9 

0.45 

208±22 

J0849-0412 

0.722 

1.54 

320±24 

1.1 

SL2S 

0.9 

0.46 

338±25 

J0849-0251 

0.274 

2.09 

276±35 

1.16 

SL2S 

0.9 

1.34 

279±35 






Table 1—Continued 


Name 

Zl 


dap [km/s] 

Oe ["] 

survey 

Oap n 

Oeff ["] 

do [km/s] 

J0850-0347 

0.337 

3.25 

290±24 

0.93 

SL2S 

0.7 

0.28 

309±26 

J0855-0147 

0.365 

3.39 

222±25 

1.03 

SL2S 

0.7 

0.69 

228±26 

J0855-0409 

0.419 

2.95 

281±22 

1.36 

SL2S 

0.7 

1.13 

283±22 

J0904-0059 

0.611 

2.36 

183±21 

1.4 

SL2S 

0.9 

2 

182±21 

J0959+0206 

0.552 

3.35 

188±22 

0.74 

SL2S 

0.9 

0.46 

199±23 

J1359+5535 

0.783 

2.77 

228±29 

1.14 

SL2S 

1 

1.13 

233±30 

J1404+5200 

0.456 

1.59 

342±20 

2.55 

SL2S 

1 

2.03 

342±20 

J1405+5243 

0.526 

3.01 

284±21 

1.51 

SL2S 

1 

0.83 

294±22 

J1406+5226 

0.716 

1.47 

253±19 

0.94 

SL2S 

1 

0.8 

262±20 

J1411+5651 

0.322 

1.42 

214±23 

0.93 

SL2S 

1 

0.85 

221±24 

J1420+5258 

0.38 

0.99 

246±23 

0.96 

SL2S 

1 

1.11 

252±24 

J1420+5630 

0.483 

3.12 

228±19 

1.4 

SL2S 

1 

1.62 

230±19 

J2203+0205 

0.4 

2.15 

213±21 

1.95 

SL2S 

1 

0.99 

219±22 

J2205+0147 

0.476 

2.53 

317±30 

1.66 

SL2S 

0.9 

0.66 

330±31 

J2213-0009 

0.338 

3.45 

165±20 

1.07 

SL2S 

1 

0.27 

179±22 

J2219-0017 

0.289 

1.02 

189±20 

0.52 

SL2S 

0.7 

1.01 

191±20 

J2220+0106 

0.232 

1.07 

127±15 

2.16 

SL2S 

1 

0.8 

132±16 

J2221+0115 

0.325 

2.35 

222±23 

1.4 

SL2S 

1 

1.12 

227±24 

J2222+0012 

0.436 

1.36 

221±22 

1.44 

SL2S 

1 

1.56 

223±22 
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Table 2. Dark energy {XCDM model and CPL parametrization) constraints obtained on the 

full 118 strong lensing (SL) sample. 


Cosmology (Sample) 


Wo 

Wl 


70 


71 

XCDMl (SL; aap) 

Wq = 

-1 45+0-54 
-^•+^-0.95 

wi = 0 

7o 

= 2.03 ± 0.06 

7i 

= 0 

XCDMl (SL; uo) 

Wo = 

1 ip:+0.56 
-‘-•■‘-^-1.20 

Wl = 0 

70 

= 2.07 ±0.07 

71 

= 0 

XCDM2 (SL; aap) 

Wo = 


Wl = 0 

70 

= 2.06 ± 0.09 

71 

= -0.09 ±0.16 

XCDM2 (SL; ao) 

Wo = 

1 qc+o.e? 

Wl = 0 

70 

9 1 Q+U.UV 

— Z.10_Q-^2 

71 

= -0.09 ±0.17 

CPLl (SL; aap) 

Wo = 


Wl = -Q.95™ 

70 

= 2.08 ±0.09 

71 

= -0.09 ±0.17 

CPLl (SL; ao) 

Wo = 

—1 on+^-54 

Wl — —1 S'^+4-85 
Wl — i.50_g 75 

70 

- 2 14+0-07 

71 

= -0.10 ±0.18 

CPL2 (SL; aap) 

Wo = 

—0 16+^-^^ 
U.iD_i 48 

Wl = -6.25+1^^ 

70 

= 2.08 

71 

= -0.09 

CPL2 (SL; ao) 

Wo = 


Wl = 

70 

= 2.14 

71 

= -0.10 

CPL2 (SN) 

Wo = 

-1.00 ±0.40 

Wl = —0.12^2.78 

□ 


□ 



^'In our fits we separately considered observed velocity dispersions aap and corrected velocity dis¬ 
persions (To, XCDMl corresponds to assumption of non-evolving power-law index 7 , while XCDM2 
assumes its evolution 7 ( 2 ;) = 70 + 7 i zi. Fixed prior of = 0.315 was assumed according to the 
Planck data. While fitting CPL parameters we assumed evolving lens mass density with 70 and 71 as 
free parameters (CPLl) and then fixed them at best-fit values (CPL2). For comparison fits of CPL 
parameters using Union2.1 supernovae data (SN) is shown. 







